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We define and explore the classical counterpart of entanglement in complete analogy with quan- 
tum mechanics. Using a basis independent measure of entropy in the classical Hilbert space of 
densities that are propagated by the Frobenius- Perron operator, we demonstrate that at short times 
the quantum and classical entropies share identical power laws and qualitative behaviors. 
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Quantum entanglement has been studied with renewed 
vigor in the recent past due to its relevance in quan- 
tum information theory p|. It has been recognized as 
an unique quantum resource that could be utilized for 
potentially powerful and qualitatively different quantum 
computers O] . It measures the degree of non-separability 
between individual parts of a composite system. That 
this non-separability may have apparently dramatic ef- 
fects was implicitly exploited in the Einstein, Podolsky, 
Rosen paradox. It appears that much of the strangeness 
of quantum phenomena follows from the existence of such 
entanglement. A classical system, composed perhaps of 
very many interacting parts, could have at all times a 
unique state assigned to these parts, in terms of their 
phase-space variables. However it is quite clear that a 
quantum state is not the counterpart of a classical phase- 
space point; although both describe the maximum pos- 
sible information we can have of the system within the 
framework of the respective theories. 

To arrive at the classical counterpart of entanglement 
we study a classical measure of non-separability of indi- 
vidual subsystems. In the spirit of entanglement we do 
not allow for transformations that may mix the states 
of the individual subsystems, i.e., the allowed transfor- 
mations must be local. A Hilbert space formulation of 
classical mechanics is the natural starting point towards 
this end. Recently, classical mechanics on the Hilbert 
space has been studied to understand certain aspects of 
quantum-classical correspondence, especially in the pres- 
ence of classical chaos |^. The principal object is the 
Frobenius-Perron (F-P) operator, that is the exponenti- 
ation of the Liouville operator, which time evolves phase- 
space densities and is the classical counterpart of the 
quantum unitary propagator. States in the Hilbert space 
are allowed phase-space densities, and the measure of 
separability is simply the separability of evolving phase- 
space densities. Thus the machinery of quantum entan- 
glement, measuring entanglement in a complex Hilbert 
space, may be applied directly within this space of clas- 
sical densities. Recent work, with the aim of bringing 
into sharper focus the unique nature of quantum entan- 
glement, has explored the classical part of correlations 
from an information theoretic point of view B . 



Another related aspect is the issue of entropy produc- 
tion, and the definition of entropy in irreversible pro- 
cesses [p| . The usual definition of entropy in an ensemble 
through the integral version of the Shannon entropy, runs 
into conceptual difficulties. As is well known it is strictly 
a constant and thus some sort of coarse-graining is needed 
to define a changing entropy. Entropies are then defined 
through partially integrated out densities, say through 
the Bogoliubov hierarchy |q|. On the other hand a von 
Neumann entropy of the classical densities, as we define 
below, will be non-negative by construction, need not be 
a constant of motion, and may be thought of as the most 
natural possible coarse-graining. A related work in this 
context which explores quantum-classical correspondence 
for open chaotic systems is in M. 

Due to the difficulty in defining a measure of entangle- 
ment for multipartite spaces, we consider a system having 
just two particles. Quantum mechanically the system is 
described by states in a bipartite Hilbert space and as is 
well known there exists a Schmidt decomposition which 
if it has more than one non-zero term implies entangle- 
ment, and the converse is also true ||l|. Now consider 
the classical mechanics as described in the phase-space 
{qi,Pi), i = 1,2. The flow is generated by a Hamiltonian 
H and the F-P operator V is defined through: 

V{t)^cM-t{ ,H}); p{q,,p,,t)=V{t)p{q„p„0). (1) 

The F-P operator is unitary. Rather than considering a 
Hilbert space of L^ functions, for simplicity we consider 
the space of L^ functions. For instance we may define the 
evolution of the square-root of a classical density. This 
also seems natural in view to our wanting to do compar- 
isons with quantum entanglement. Let {-F'm.n('Z,p)} be 
an orthonormal basis set in such a Hilbert space, which 
we will write compactly as {\m, n)}. We consider the case 
where there is a denumerably infinite number of such ba- 
sis functions labeled by the pair {m,n). A pair is used 
to facilitate defining a basis in each Hilbert space corre- 
sponding to each phase-space variable. Thus we can con- 
sider the classical Hilbert space as being an outer product 
of four Hilbert spaces. However in the following we will 
consider it as a product space of two Hilbert spaces, each 
corresponding to one degree of freedom. 



Expanding any density p in the four-dimensional phase 
space in terms of these functions we may write: 

\p) = ^a(mi,ni;m2,n2)|mi,ni) |m2,n2), (2) 

where the a(; ) are expansion coefficients. Now the "re- 
duced density matrix" is simply the partial trace: 



p\ = Tt2{\p){p\), 



(3) 



while the von Neumann entropy of this reduced density 
matrix is the entropy of classical entanglement and is 
given by: 



Sc^Tr,{p\\ogip\)). 



(4) 



The analogy to the quantum definitions are evident and 
the existence of a Schmidt decomposition follows. The 
entropy Sc can be zero or positive, if it is zero there is no 
entanglement, and the densities are separable into each 
degree of freedom, while if it is positive this is not pos- 
sible. The basis independence of the entropy Sc follows 
from the two trace operations, and thus it represents a 
physically meaningful quantity, which by its construction 
ought to qualify as the classical counterpart of quantum 
entanglement. 

We now provide evidence for this by means of an ex- 
ample. We choose two coupled standard maps for the 
possibility of studying potentially interesting connections 
between deterministic chaos and entanglement |8|. In 
fact the motivation for the present work is derived from 
a previous work g| where we studied a similar system 
which showed that classical chaos enhances quantum en- 
tanglement. Furthermore we showed that some Schmidt 
vectors were "scarred" by projections of certain classi- 
cal periodic orbits. A study of coupled kicked tops has 
also shown that quantum entanglement increases linearly 
in time at a rate which is the sum of the classical posi- 
tive Lyapunov exponents ||l^ . These are indications that 
quantum entanglement is not entirely immune to the lure 
of the classical. 

Consider the classical map defined on the four-torus 
r2 X T^: 



<l'i = (li+Pi (modi) 
p[^Pi-dV/dq[ (modi). 



(5) 



where i = 1, 2 and the potential V is 
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■ cos(27rgi) cos(27r(72) (6) 



This is a symplectic transformation on the torus T** and 
may be derived from a kicked Hamiltonian in the stan- 
dard manner S . This is somewhat different from the pre- 
viously studied four-dimensional generalizations of the 



standard map, and has been chosen to maximize symme- 
tries. The parameter b controls the interaction between 
the two standard maps, while Ki determine the degree of 
chaos in the uncoupled limit. The advantage in studying 
such maps is our ability to write explicitly the relevant 
unitary operators, both classical and quantal. 

The classical unitary F-P operator over one iteration 
of the map is 



V^Vl®V2Vb: 



(7) 



where Vi is the F-P operator for the standard map on T^ 
and acts on the single-particle Hilbert spaces, while Vb is 
the interaction operator on the entire space. We use as 
a single-particle basis the Fourier decomposition: 

Fm^n[q,p) ^e:-x:p{2m{mq + np)) . (8) 

Then the matrix elements of the F-P operators are: 

{'m'^,n[\'P^\mi,n.i) = im[--m, (^^"D '5n,-m,,n;, (9) 
and 



H+ 



m[,n[;m'2,n'2\Vb\{mi,ni;m2,n2) = 

Ji- ( "'""„ ) Sni.n',Sn2.n'.,, (10) 



6(ni+7i2)\ (b{ni—n2) 



where 



l± 



((mi 



l) ± (ma - mf,)) , 



and the delta functions are Kronecker deltas. The orders 
of the Bessel J functions are restricted to the integers. 

If the initial density is described by a vector (Eq. (0)) 
|p(0)) then 



\p{T))^V^\pm 



ill) 



is the density after a time T. The same may of course 
be equivalently obtained by means of the Liouville equa- 
tion. The number of classical components excited in- 
creases with time in general including increasingly larger 
frequencies. This is indicative of the fine structure that is 
being created in phase-space by the stretching and fold- 
ing mechanism. The frequencies involved appears to in- 
crease exponentially in time for chaotic systems, while 
we may expect a polynomial growth for integrable or 
near-integrable systems. The Arnold cat map provides 
a completely chaotic exactly solvable model where this 
exponential increase is readily seen |ll|, while the free 
rotor, corresponding to the case Ki ~ for the standard 
map on T^ is an exactly solvable integrable model where 
a linear growth is apparent. With the flow of probabil- 
ity moving inexorably outward, for chaotic systems, the 
infinite dimensional Hilbert space of classical mechanics 
provides a setting for an approach to equilibrium, with 



every smooth density tending towards the invariant den- 
sity in a weak sense Q . 

The quantization of the symplectic transformation in 
Eq. (^ is a finite unitary matrix on a product Hilbert 
space of dimensionahty iV^, and N = l/h, where h is a 
scaled Planck constant. The classical limit is the large 
N limit. The quantization is straightforward as there 
exists a kicked Hamiltonian generating the classical map 
|g[ . The quantum standard map on the two-torus in the 
position representation is 



/ /1,-ri X exp(— i7r/4) f.^rKi ,2tt . 1,, 

{n',\U,\m) = \^^ exp KiV^cos( — (n, + -)) 



X exp( j^{n,-n[f 



(12) 



The position kets are labeled by n = 0, N— 1 and the po- 
sition eigenvalues are (n -I- l/2)/N while the momentum 
eigenvalues are m/N, m = 0, . . . , N — 1. The quantum 
phases have been chosen to maximize symmetries, and N 
is taken to be even. The four-dimensional quantum map 
is but a simple extension: 



U = Ui®U2 Ub, 



(13) 



which is the quantum equivalent of the classical F-P op- 
erator in Eq. (0) and the interaction operator is 



{n[,n'2\Ub\ni,n2)^ exp 



-iN — cos — [n-i -\ — 



C0S(— (^2 + -) 



^rii ,n' ^n2 ,n' 



(14) 



J7 is a unitary matrix and will induce mixing between the 
two subsystems. Some properties of such unitary matri- 
ces especially in connection with quantum entanglement 
has been studied previously . 

For the initial state a product state composed of coher- 
ent states is a natural choice, as the corresponding initial 
classical density becomes apparent. We choose 



|V(0)) = |00) ® |00), 



(15) 



where |00) is a coherent state that is peaked at the clas- 
sical point {q = 0,p = 0) Q, which also happens to be 
a classical fixed point for the dynamics. The initial clas- 
sical density will be a circular (periodicised) gaussian of 
width a around the origin. We choose 

ao{m, n) = (87ro-2)i/2 exp{-4Tr^a^{m^ + n^)), (16) 

to be the Fourier components of such a (approximate) 
periodicised gaussian. The initial density's components 
are then specified by the outer product 

ao(mi,ni;m2,n2) = ao(mi, ni) ao(m2, 712). (17) 



For a ^ l/vN the Husimi distribution of the quantum 
state is similar to the (square of the) classical initial den- 
sity. For quantum-classical correspondence we require 
large N or small a. 



We time evolve the initial states, both quantum and 
classical, using the respective propagators and calculate 
the entanglement in each. Due to the infinite dimension- 
ality of the classical Hilbert space and the at least poly- 
nomially increasing frequency components, the classical 
calculations in a truncated Hilbert space lead to rapid 
loss of accuracy. Thus this work, whose main intention 
is to compare the classical entanglement as defined in 
Eq. (0) with quantum entanglement, merely touches the 
surface of interesting results. We put Ki = K2 = and 
turn on the interaction parameter b only slightly. Thus 
the uncoupled systems are free rotators, while the cou- 
pled system is near-integrable. For sufficiently large val- 
ues of the interaction, full fledged chaos seems to develop, 
a case which we do not discuss further here. 



First we look at entanglement produced for very short 
times, as this is numerically easily accessible. We study 
this entanglement as a function of the interaction. In 
Fig. n^ is shown entanglement produced at times 1 and 2. 
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FIG. 1. Classical and quantum entanglement after time 
T = 1 (top) and T = 2 (bottom) as a function of the interac- 
tion for the initial state described in the text. The quantum 
calculations correspond to A*' = 50, while the classical initial 
state has cr = 0.1. 



We find approximate power laws that are nearly the same 
for both classical and quantum entanglement. The best 
correspondence is naturally at the shortest time, T = 1, 
when a power law of the form S ~ b^'^ holds in both 
cases, while at time 2 the exponents are not so close. 
However at this time the power law seems to be approx- 
imate and the deviations from the straight line is seen to 
have common features in both the classical and quantum 
cases. The N or h dependence of the quantum entan- 
glement does not seem to be a power law at these short 
times and we do not pursue this here, thus we do not 
compare the absolute values of the two entanglement. 

In Fig. g we plot both the entanglements as a function 
of time for a fixed value of the interaction. We see qual- 
itatively similar behaviour. The quantum entanglement 
tend to saturate at higher values for larger N, while the 
classical entanglement seems to be continuously increas- 
ing. This is however difficult to see numerically due to 
increase of errors, for instance at the highest time shown 
in the figure, T = 6, the normalization deteriorated from 
unity to about .9948. More elaborate work is needed to 
establish the properties of classical entanglement. 
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Time T 
FIG. 3. Quantum entanglement as a function of time for 

the same initial state as in Fig. hi. Ki = 6.0, K2 — 5.0, while 

The interaction is fo = 0.001. The linear regime is seen to 

scale as A^^. 



The linear entropy increase indicates an exponentially 
rapid increase in the number of single-particle basis states 
used. Thus we conjecture that this regime is in fact a 
classical one, with this being a reflection of the expo- 
nential increase in the number of basis states explored 
in the classical Hilbert space. The rate being propor- 
tional to the sum of positive classical Lyapunov expo- 
nents seems to suggest that the classical entanglement 
production rate is proportional to the classical K-S en- 
tropy. The subsequent saturation is a quantum effect, 
while we may expect the classical entropy to continue to 
linearly increase. We further conjecture that the differ- 
ence between chaotic and near-integrable or integrable 
systems will be reflected in the classical entanglement 
SciT) behaving at large times as either T or log(T), re- 
spectively. Much larger classical calculations will be able 
to prove or disprove these statements. 



Time T 
FIG. 2. Classical and quantum entanglement as a function 

of time for the same initial state as in Fig. Fjl The interaction 

6 = 0.05. 



We focus finally on an interesting regime where both 
the subsystems are chaotic, while the interaction is weak. 
In this case a linear quantum entropy increase was found, 
with the rate being proportional to the sum of positive 
Lyapunov exponents RJlQl- ^^ ^^'^^ this linear regime is 
time-bound and gives way to saturation, while the en- 
tropy in this linear regime scales as N'^ (^2). This is 
illustrated in Fig. for the system we consider in this 
paper. 
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